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Abstract 

We analyze the stability of soliton solutions in a Cliern-Simons-CP(l) model. We 
show a condition for which the soliton solutions are stable. Finally we verified this 
result numerically. 
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1 Introduction 

The two dimensional CP{n) sigma model was introduced in the late seventies [UlSlll], 
in the search of understanding the strong coupling effects in QCD. This model captures 
several interesting properties, many of them present in four dimensional QC-D [ll[5l|6l[7|. 
Whereas in four dimensional QCD is difficult to demonstrate the existence of these 
properties, in two dimensional CP{n) sigma model it becomes comparatively simple. 
An important issue related to this type of models concern to the existence of soliton 
type solutions. For the simplest CP{1) model topological solutions have been shown 
to exist[8]. Nevertheless, the solutions are of arbitrary size due to scale invariance. As 
argued originally by Dzyaloshinsky, Polyakov and Wiegmann[9] a Chern-Simons term 
can naturally arise in this type of models and the presence of a dimensional parameter 
could play some role stabilizing the soliton solutions. A first detailed consideration 
of this problem was done in Ref.[lOj where a perturbative analysis around the scale 
invariant solutions (i.e no Chern-Simons coupling k = 0) showed that the solutions 
were pushed to infinite size. More recently, in Ref.[llJ, a nonperturbative analysis 
of the solutions was done, showing that the Chern-Simons-CP(l) system, without a 
potential term, admits only trivial solutions in R^. 
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This paper pretends to be a continuation of the work and a generahzation 
of the results obtained there. We will show that a Chern-Simons-CP(l) model with 
a potential term, which was proposed in the reference [12], presents a stable soliton 
solution if there is a critical radius Sc such that the equality 

/ (fx = [ (fx V (1) 

Jds Jds 

is satisfied for all radius S > where here the subindex Ds indicates that the region 
of integration is a disc of radius S and the letters V and B represent the potential term 
and the magnetic field. 



2 The model 

We begin by considering a (2 + l)-dimensional Chern-Simons model coupled to a com- 
plex two component field n{x) described by the action 



D 

The subindex D indicates that the region of integration is a disc D of radius i?[TT]. 
Here = 9^ — iA^ (/i = 0, 1, 2) is the covariant derivative and V is the potential term 
to be determined later. The term Scs is the Chern-Simons action given by 



S = Scs+ d'^xlDf.nl' + V (2) 



Sc 

where 



K [ (fxe^^'^PA^d.Ap (3) 

JD 



F^u = d^A, - d,A^ (4) 

The metric signature is (1, —1, —1) and the two component field n{x) is subject to the 
constraint n)n = 1. The constraint can be introduced in the variational process with 
a Lagrange multiplier. Then, we extremise the following action 

S = Scs + j d^xi^lD^nl"^ + V + (5) 
The variation of this action yields the field equations 

dV 

D.D^n - — + An = (6) 
t^e^upF'^P = -J^ = i[v)D^n - n{D^n)^] (7) 



From the first of these equations we get A = —n^ (D^D^^n — ), so that 



(9F / + dV \ 



The time component of Eq. ([7]) 

2kFi2 = - Jo (9) 

is Gauss's law of Chern-Simons dynamics. Integrating it over the entire plane one 
obtains the important consequence that any object with charge Q = cPxJq also 
carries magnetic flux ^ = B(Px [131 EH ES]: 

^ = -^Q, (10) 

where in the expression of magnetic flux we renamed F12 as B. 
Defining the stress-tensor as|12j 

V = iD^n)^D,n + {D,n)^D^n - g^, (^{D^n)^D"n - , (11) 

we get 

E = j d'^xi^K'^B'^ + \Din\'^ + , i = l,2, (12) 

which is the expression of the energy functional for the static field configuration. 

Here, It is convenient to specified the potential term V. Following the reference 
[I2] we define the potential as 



V{n) =r](l-n^{a3n)^ (13) 

Where 

is the third Pauli spin-matrix and r/ is the coupling strength . 

Let us consider, now, the following ansatz with cylindrical symmetry for the A'^ 
soliton solutions: 

/ cosf^le^^-^ \ 

n{cl>, r) = i 2n{'-^) ) ' ^"^^"^ = ""^"^ ' = ° ' ^^^^ 

In terms of this ansatz the energy (jl2p reads as 

cos'(^) + a^r) + r?(l - cos(0(r)))) , (16) 



'N^ 2Na{r)\ _2/^('') 



^2 



3 



whereas the field equations become 



2 dra{r) a{r) a{r) _ 2/^(^)\ ^ 



rdr{rdre{r)) + (^iV^ _^ 2iVra(r)) sin(6'(r)) = r?2r2 sm(6'(r)) , (18) 
In order to ensure the regularity of the field at the origin, we impose 



lim e(r) = TT , lim a(r) = (19) 

On the other hand, the conditions at the boundary of the disk are in principle more 
general. This is because the length of the radius R is arbitrary. However, if the size 
of the disc becomes infinite, then, we must impose boundary conditions ensuring finite 
energy, that is 

lim e{r) = 0, lim a{r) = (20) 

r— >oo r— >-oo T 

Therefore, it is convenient to use this boundary condition independently of the length 
of the radius. So, we fix the boundary conditions at R to be 

lim 9(r) = , lim a(r) = (21) 

r^R r^R r 

this boundary conditions imply the quantization of the magnetic fiux 

<^ = 2n rrdr^^^^^^ = -27rN (22) 
Jo r 

If the solutions of (fTTj) and ([TS]) exist their scale must be set by the quantity k. Following 
Ref.|16j. we introduce the dimensionless quantities 

A = Ka, s = - (23) 

K 

in terms of which (flTl) and (fTSi) become 



dlA^^ ^_^ = cos2-- 24 



sdsisdsO) + (n^ + 2NsA^ sin(0) = 7j2s'^ sm{e) (25) 
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The energy functional (jl6p in terms of these new variables reads as 



E{S) = 2Tr j\ds(^(^^ + dsA^ +^{dsef 

+ + ^) cos^(^) + A' + Vil - cos(0))) (26) 
For the origin we choose the following boundary conditions, 

hm 6* = TT , lim ^ = (27) 

while for the boundary 5* = R/k we choose, 

1- . ^ 
lim A = — — 

s^S S 

lim 6* = (28) 

s^S 

3 Stability analysis 

In this section we analyze the stability of soliton solutions corresponding to the equa- 
tions dM]) and ([25]). 

Consider the following configuration defined in the interval of length AS* 

A,sis) = ^^, 0^sis) = es{j) (29) 



Here A is a positive real number such that A > 1 and the configurations (j29p satisfy 
the boundary conditions 



lim Oxsis) = vr , lim Axs{s) = (30) 

s->-0 s— !>0 



-AT 

lim Oxsis) = , hm Axs{s) = — (31) 

s-s>A5 s^XS Ad 

We can evaluate the energy functional (|26p for the configuration (j29p in an interval of 
length AS* 

E{XS) = 27T £\ds(^ l^^^^ + dsAxsis)^ +\{dsexsf 
I ^ s J ^^^^> + ~ cos(^As)) j (32) 



We denote the solution corresponding to the interval AS" as Axs{s) and 6xs{s) and its 
energy as E{\S). Since the configuration (I29p satisfy the same boundary condition as 
^\s{s) and 9xsis), we have that 

E{XS) < E{\S) (33) 
Under the transformation s = xX the functional ([32]) becomes 

E{XS) = 2^ j\dxQ^ + dMx)j + \{dJs? 

Since there are evaluated in the same interval, we can compare this expression with the 
formula (j26p . For this purpose we can look for the values of A for which the equality 
E{XS) = E{S) is held. Using the equations ([26]) and (l34l) . we obtain 

= xdx ( (^^^ + dMx)) ' + - cos(0s))) (35) 

Renamed 



X^ = oj 

s 

xdxr]{l — cos{6s)) = a 



\dx ( + d^Asix)] = b, (36) 



/o \ X 
the expression ([35]) reduce to 



Lo'^a -{a + b)uj + b = (37) 
The roots of this polynomial are 

a; = 1 , a; = - (38) 
a 



and then the possible values of A are 



A = l, Ae = A/- (39) 
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Now, suppose that b > a and choose the values of A satisfying 1 < A < y ^. It is not 
difficult to see that for this values of A the following inequality is held 



< I' xdx ( (^^^ + d.Asix)^ ' + r?(l - cos(0s))) (40) 
and therefore 

E{\S) < E{S) (41) 

Comparing (I33p and (I4ip we have that 

E{XS) < E{S) (42) 

that is, the energy decreases when we enlarge the interval S. This process is valid only 
for 1 < A < ^J^■ We can repeat the proses by considering an interval of length XcS 
instead of the interval S. Now, we have 



\cS 

xdxrj{l - cos{6\^s)) = «i 



^^\dx{^^^ + d,A^^s{x)^^ = h, (43) 

again, if oi < 61, we can conclude that the energy decreases when the interval is 
enlarged. That is 



EiXiX^S) < E{XcS) (44) 

where 1 < Ai < The process can be repeated for successive intervals provided 

that b > a in each of this intervals. Therefore if the relation 



J^^ xdx7]{l - cos{e)) < xdx [^^^ + d^As{x)^ (45) 

is held for all interval, the energy decreases indeterminably as 5 ^ 00, and thus there 
are no finite size soliton solution in R^. 
Now, suppose 

rpSi rpSi / A Q (x) \ ^ 

xdxvil - cosiBps,)) > xdx iy^^^ + d,Aps,ix)j (46) 



where indicates the length of an arbitrary interval and /? is a real number such that 
p > 1. Consider the following configuration defined in the interval of length Si 



As, (s) = pAps, (ps) , es, is) = Ops, [ps) (47) 

where Aps,{s) and 9pSi{s) are the solutions of the field equations in the interval pSi. 
In virtue of equation ([27|) and ([28]) . Aps,{s) and Ops,{s) must satisfied 



hm Ops, (s) = vr , Hm Aps, (s) = (48) 

s— >-0 s— >-0 



hm Ops, (s) = , hm Aps, [s) = — (49) 

s— s-^pbi pji 

Then, the configuration ([^7]) is subject to the following boundary conditions 



lim 03. (s) = TT , Hm As, (s) = , (50) 

s->-0 s^O 



lim OsAs) = 

hm ^5i(s) = -P^ = ^ (51) 
s~^Si pbi bi 

The solutions of the field equations in the interval Si, which we denote by Asj^{s) and 
9si{s), also satisfied the boundary conditions (fSOl) and (fST]) . Therefore, we have 

^(5i)<^(5i) (52) 
where E{Si) is the energy corresponding to the solution ^45^(5) and 9si{s), and 

E{Si) = 27r l\ds(^(^^^ + dsAsAs?j +\{dsOs,f 

+ (5 + ^^^)cos^(%) + i|.(^)+.(l-cos(^~.^ (53) 

Using the configuration (H7|) and under the transformation x = ps, the expression ([53]) 
reads as 



E(5i) = 27r iyxdx[p'(^^^^+d.ApsAx)^ +\{d.^~'' 



'N^ 2NApsAx 

11 o ~r 



cos2(^) + AlsA^) + ^(1 - cos(0,5j)) (54) 
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The energy for the solutions in the interval pSi is 



^) cos2(^) + Als, (x) + vil - cos(0p5j)) (55) 



2NAr,s. (x) 



+ 1- + ^^^ 

X 

Since ()54p and (j55p are evaluated in the same interval we can compare this formulas. 
In this case we look for the values of p for which the equality E{Si) = E{pSi) is held. 
Following the same steps that we do previously we arrive to the polynomial 

uj'^a - (a + 6)a; + 6 = (56) 

However, in this case a and h are different from the constants present in the polynomial 
(1371). In fact we have 



p^ = OJ 



xdx \-d:,Aps^{x^ 







X 
pSi 



xdxri{l — cos{9pSi)) = b (57) 
As in formula (j39|) the roots of the polynomial produce the following values of p 







P = l, Pc = \- (58) 
V a 

Again, it is easy to show that choosing the values of p satisfying the relation 1 < p < 
- we have 



EiSi) < E{pSi) (59) 

and therefore in virtue of (1520 

E{S{) < E{pSi) (60) 

We can repeat the process by choosing pcSi instead of Si. In that case the formulas 
HD and dll]) read as 

xdxr,{l-cos{ep,p^s,))> xdxy ''"^^'^ ' + d^Ap,p^sA^)\ (61) 
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^Pi5i (s) = /5i^pipc5i ipis) , 9p^si (s) = Op^p^si ipis) (62) 
and finally we arrive to 

E{p,Si) < E{pip,Si) (63) 

Of course the process can be repeated indefinitely provided that 6 > a in all intervals. 
This implies that if the relation 

j^^ xdxr]{l - cos(05i)) > jj' xdx [^^^ + d^As^{x)^ (64) 
is valid for all interval ^i, the energy increases. Certainly, we have 

< E{pSi) - E{Si) < E{pSi) - E{Si) , (65) 



where 



E{pSi)-E{Si) = {l-p^)a+{l~\)b (66) 
The roots of this equation are 

p = l, P=\J\ (67) 

Since 1 < p < \J~^ for all interval, the energy increases indeterminably as 5 — )• oo. The 

stability only can take place if \E[pSi) — E{Si) \ and this implies that the interval 

1 < p < ^J~^ must be contracted to a point, that means a = b and then p = 1. 

We conclude that if a soliton solution exist, then there is a critical radius Sc such 
that the equality 



J xdxr]{l — cos{6^)) = J xdx y~ ^ dxAg{x)j 

is verified for all radius S > Sr- 



4 Numerical Results 

The field equations ()24p and ()25p can be resolved numerically. From the numerical 
point of view one solve the field equations in a finite disc and then the analysis of the 
solution in an infinite plane can be carried out by increasing successively the radius 
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Figure 1: The energy density, without a potential term, as a function of scaled radial 
coordinate s for different disc sizes, from top to bottom, S—30, 50, 60, 70, 100, 150. 
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Figure 2: The energy density, including a potential term, for different disc sizes, from top 
to bottom, 5=10, 30, 40, 50, 60. The solution was calculated using a coupling constant 
Tj = 10"^ 
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of disc. In order to solve numerically the field equations for arbitrary disc sizes, the 
relaxation method is suitable. In Fig|T]we show the behavior of the energy density, for 
N = 1 and S = 30, 50, 60, 70, 100, in the especial case in which there is no potential 
term and therefore the inequality ()45p is held for arbitrary disc sizes. This fact is 
due to the constant value of the magnetic flux, which implies B ^ and therefore 
f^^(PxB'^ > for all finite disc. This corresponds to problem analyzed in Ref.|llj. 
The figure clearly show that as the disc size increases the energy density tends to zero 
and then the solution is not stable in a finite disc. These results are in concordance with 
our previous analysis which point out that if the inequality (|45p is held for arbitrary 
disc sizes, then the energy tends to zero. 

In Figi2]show the energy density, for = 1 and for disc sizes S = 10, 30, 40, 50, 60, 
with the inclusion of a potential term. We can see that the energy density is highly 
instable for small disc sizes the size (i.e. S = 10, 20) and becomes more stable as S 
increases. In fact the energy density corresponding to = 50 and S = 60 are practically 
indistinguishable in the graph and becomes completely stabilized for S = 60. From 
theoretical point of view we have that the inequality (I45p is held for < 60 and it 
becomes an equality for S > 60. So, Sc = 60 is the critical radius from which the 
equality ([681) is verified. 

5 Conclusion 

In summary we have studied the classical solution of the Chern-Simons-CP(l) model 
with a potential term. Specifically we have shown that if the soliton solution exits, 
then there is a disc Ds^ such that 



is satisfied for all disc Ds > Ds,.- In addition we resolved numerically two situations. 
In the first situation we analyzed the model without a potential term, showing that 
the energy decreases as S ^ oo , which implies the instability of the solutions. As 
a second case we analyzed the model with a potential term. We showed the energy 
decreases for S < 60 and becomes stable for S > 60. 
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